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Abstract 

In  this  paper,  we  develop  nonlinear  constitutive  equa¬ 
tions  and  resulting  system  models  quantifying  the  non¬ 
linear  and  hysteretic  field-displacement  relations  in¬ 
herent  to  piezoceramic  elements  employed  in  present 
nanopositioner  designs.  We  focus  specifically  on  piezo¬ 
ceramic  rods  utilizing  ^33  motion  and  piezoceramic 
shells  driven  in  d^\  regimes,  but  the  modeling  frame¬ 
work  is  sufficiently  general  to  accommodate  a  variety  of 
drive  geometries.  In  the  first  step  of  the  model  develop¬ 
ment,  lattice-level  energy  relations  are  combined  with 
stochastic  homogenization  techniques  to  construct  non¬ 
linear  constitutive  relations  which  accommodate  the 
hysteresis  inherent  to  piezoceramic  compounds.  Sec¬ 
ondly,  these  constitutive  relations  are  employed  in  clas¬ 
sical  rod  and  shell  relations  to  construct  system  models 
appropriate  for  presently  employed  nanopositioner  de¬ 
signs. 

1.  Introduction 

Pieozoceramic  materials  have  been  employed  as  ac¬ 
tuators  in  scanning  tunneling  microscopes  (STM)  and 
atomic  force  microscopes  (AFM)  since  their  inception 
due  to  their  high  set  point  accuracy,  large  dynamic 
range,  and  relatively  small  temperature  sensitivity  [8]. 
More  recent  investigations  have  focused  on  the  design 
of  AFM  stages  for  employment  in  applications  ranging 
from  nanoconstruction  to  the  development  of  nuclear 
magnetic  resonance  microscopes  (NMRM)  with  the 
goal  of  detecting  single  electron  spins  [6,  13,  20].  How¬ 
ever,  a  fundamental  challenge  when  employing  piezoce¬ 
ramic  actuators  even  at  the  low  drive  levels  required 
for  nanopositioning  is  the  quantification  and  accom¬ 
modation  of  hysteresis  and  constitutive  nonlinearities 
inherent  to  the  materials  as  depicted  in  Figure  1. 

For  certain  drive  regimes,  the  hysteresis  and  constitu¬ 
tive  nonlinearities  can  be  mitigated  through  either  the 
drive  electronics  or  feedback  loops  incorporated  in  the 
software.  As  detailed  in  [10,  11],  the  use  of  charge  or 
current  controlled  amplifiers  can  essentially  eliminate 
hysteresis.  However,  this  mode  of  operation  can  be 
prohibitively  expensive  when  compared  with  the  more 
commonly  employed  voltage  controlled  amplifiers,  and 
current  control  is  ineffective  if  maintaining  DC  offsets 
as  is  the  case  when  the  x-stage  of  an  AFM  is  held 
in  a  fixed  position  while  a  sweep  is  performed  with 


the  y-stage.  For  low  scan  rates,  PID  or  robust  con¬ 
trol  designs  can  be  employed  to  accommodate  hystere¬ 
sis  [4,  14].  However,  at  the  high  scan  rates  required 
for  real-time  product  diagnostics  or  monitoring  of  bi¬ 
ological  processes,  increasing  noise-to-data  ratios  and 
diminishing  high-pass  characteristics  of  control  filters 
preclude  a  sole  reliance  on  feedback  laws  to  eliminate 
hysteresis.  This  motivates  the  development  of  control 
designs  which  incorporate  and  approximately  compen¬ 
sate  for  hysteresis  through  model  inverses  employed  ei¬ 
ther  in  feedback  or  feedforward  loops. 

In  this  paper,  we  develop  models  for  this  purpose 
through  a  two  step  process.  In  the  first  step,  Gibbs 
energy  relations  at  the  lattice  level  are  minimized  to 
provide  models  for  the  local  average  polarization  gen¬ 
erated  by  an  applied  field.  The  effects  of  material  non¬ 
homogeneities,  polycrystallinity,  and  variable  effective 
fields  are  subsequently  incorporated  through  stochas¬ 
tic  homogenization  techniques  to  provide  macroscopic 
constitutive  relations  which  quantify  the  nonlinear  and 
hysteretic  field-polarization  and  field-strain  behavior  of 
the  materials.  In  the  second  step  of  the  development, 
these  constitutive  relations  are  incorporated  in  classical 
rod  and  shell  theory  to  provide  system  models  which 
quantify  the  displacements  generated  by  prototypical 
nanopositioning  designs. 

To  place  this  modeling  framework  in  context,  we 
compare  it  with  the  domain  wall  model  employed  in 


Figure  1.  Quasistatic  relation  between  the  input  field 
E  and  displacements  generated  by  a  PZT  stacked  ac¬ 
tuator  in  an  AFM  stage. 
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[15]  and  Preisach  model  developed  in  [16]  for  charac¬ 
terizing  the  hysteretic  behavior  of  nanopositioners.  For 
certain  choices  of  the  Gibbs  energy,  the  local  polariza¬ 
tion  model  in  the  present  theory  is  the  same  as  the 
anhysteretic  relation  in  the  domain  wall  theory.  At  the 
macroscopic  level,  however,  the  two  models  differ  sig¬ 
nificantly  -  from  the  perspective  of  implementation,  the 
present  theory  guarantees  the  closure  of  biased  minor 
loops  whereas  this  is  accomplished  in  the  domain  wall 
model  only  if  a  priori  knowledge  of  turning  points  is 
available.  Furthermore,  it  is  illustrated  in  [18]  that  the 
present  theory  provides  an  energy  basis  for  extended 
Preisach  models  which  enforces  deletion  but  not  con¬ 
gruency.  Through  the  energy  derivation,  the  proposed 
theory  inherits  several  of  the  advantages  associated 
with  extended  Preisach  frameworks  while  avoiding  the 
disadvantages  associated  with  temperature,  frequency 
and  stress-dependent  parameters. 

We  consider  two  motivating  geometries  for  the  model 
development.  The  first  encompasses  PZT  rods  of  the 
type  employed  in  AFM  stages  having  the  design  de¬ 
picted  in  Figure  2(a).  In  this  case,  strains  are  generated 
through  ^33  mechanisms  which  translate  to  longitudi¬ 
nal  rod  displacements.  Appropriate  1-D  constitutive 
relations  are  developed  in  Section  2.1  and  the  corre¬ 
sponding  system  model  is  constructed  in  Section  3.1. 
The  second  geometry  is  comprised  of  a  PZT  shell  in 
which  longitudinal,  or  z,  displacements  are  generated 
through  ^31  mechanisms  -  see  Figure  2(b).  The  2-D 
constitutive  relations  and  system  model  for  this  case 
are  respectively  developed  in  Sections  2.2  and  3.2. 

2.  Constitutive  Relations 

When  modeling  the  constitutive  behavior  of  piezo¬ 
ceramic  rods  and  shells,  we  assume  linear  stress-strain 
and  strain-displacement  relations  in  accordance  with 
classical  theory.  Furthermore,  both  classical  theory, 
[3,  9],  and  recent  experiments  have  demonstrated  that 
in  the  biased  regimes  used  to  obtain  bidirectional 
strains,  the  relation  between  the  polarization  P  and 
strains  e  is  approximately  linear,  and  we  retain  that 
assumption  throughout  the  development.  At  all  drive 
levels,  however,  the  map  between  applied  voltages  V  or 
fields  E  and  the  polarization  is  nonlinear  and  hysteretic, 
and  it  is  the  quantification  of  these  inherent  properties 
which  differentiates  the  proposed  models  from  classical 
linear  theory  for  piezoceramic  rods  and  shells. 

2.1.  1-D  Constitutive  Relations 

The  1-D  constitutive  relations  necessary  for  charac¬ 
terizing  the  strains  generated  by  PZT  rods  follows  from 
the  general  theory  developed  in  [19]  for  quantifying 
the  hysteresis  and  constitutive  nonlinearities  inherent 
to  ferroelectric  compounds.  We  summarize  here  those 
aspects  of  the  theory  pertinent  to  the  development  of 
PZT  stages  employed  in  nanopositioning  devices. 


x-y  actuator 


z  actuator 


(b) 

Figure  2.  Actuator  configurations  employed  for  sam¬ 
ple  positioning  in  AFM.  (a)  Stacked  actuators  em¬ 
ployed  as  x-  and  y-stages,  and  (b)  cylindrical  PZT 
transducer. 

It  is  illustrated  in  [19]  that  for  stress-free,  fixed  tem¬ 
perature  conditions,  the  piecewise  quadratic  Helmholtz 
energy  relation 

r  \v{p+pr)2  ,p<~Pi 

il>(P)  =  |  -  PR )2  ,  P>  Pi  (l) 

[  Hpi  -  Pr)  (w  -  Pr)  >  lpl  <  Pi 

quantifies  the  internal  energy  derived  under  the  as¬ 
sumption  that  dipoles  are  either  aligned  with  the  ap¬ 
plied  field  or  diametrically  opposite  to  it.  As  depicted 
in  Figure  3,  Pj  and  Pr  respectively  denote  the  inflec¬ 
tion  point  and  polarization  at  which  the  positive  local 
minimum  of  ip  occurs.  The  point  Pr  is  also  the  local 
remanence  polarization  at  the  domain  level.  Finally, 
the  fact  that  rj  is  the  reciprocal  of  the  slope  in  the  E- 
P  relation  after  switching  can  be  utilized  to  determine 
initial  parameter  values  when  constructing  the  model 
for  a  given  piezoceramic  compound  and  application. 

To  incorporate  elastic  interactions  and  ferroelectric 
coupling,  we  also  consider  the  extended  Helmholtz  re- 
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Figure  3.  (a)  Helmholtz  energy  V’  and  Gibbs  energy  G 
for  a  =  0  and  increasing  fields  E.  (b)  Dependence  of 
the  local  polarization  P  on  the  field  E  at  the  lattice 
level  in  the  absence  of  thermal  activation. 


lation 

</>e(P,  e)  =  ^(P)  +  \Ype2  -  Ym1£P  (2) 

where  Yp  denotes  the  Young’s  modulus  at  constant 
polarization  and  7  is  an  electromechanical  coupling  co¬ 
efficient. 

To  accommodate  applied  stresses  and  fields,  we  em¬ 
ploy  the  Gibbs  energy  relation 

G(E,  P,  e)  =  +  \yP£2  -  YM"/eP  -EP-ae  (3) 


which  incorporates  the  electrostatic,  electromechanical 
and  elastic  energy  relations. 

For  operating  regimes  in  which  thermal  relaxation 
mechanisms  are  significant,  the  local  average  polariza¬ 
tion  at  the  lattice  level  is  quantified  by  balancing  the 
thermal  and  Gibbs  energies  through  the  Boltzmann  dis¬ 
tribution 

H(G)  =  Ce~GV/kT  (4) 

where  /1  specifies  the  probability  of  achieving  an  energy 
level  G,  C  is  a  constant  of  integration,  V  is  a  reference 
volume  and  k  is  Boltzmann’s  constant.  For  brevity, 
we  omit  this  case  here  and  refer  the  reader  to  [19]  for 
details  concerning  the  extension  of  the  constitutive  re¬ 
lations  to  include  thermal  relaxation  mechanisms. 

For  regimes  in  which  thermal  relaxation  is  negligible, 
the  local  average  polarization  P  at  the  lattice  level  is 
determined  by  the  necessary  conditions 
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This  yields  the  piecewise  linear  relation 
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(5) 


(6) 


for  the  respective  cases  {r(t)  =  0},  {r(f)  ^  0  and 
E(maxr(t))  =  —Ec},  { r(t )  ^  0  and  £(maxr(f))  = 
Ec}.  The  local  coercive  field  is  given  by 

Ec  =  v(Pr  ~  -Pr)-  (7) 

The  transition  points  r  are  specified  by 


T{t)  =  {t€  (0,  Tf\  I  E(t)  =  -Ec  or  E{t)  =  Ec} , 


and 


[P(E,e-,Ec,Om  = 


denotes  the  initial  dipole  orientation  for  respective  ini¬ 
tial  fields  {P(0)  <  —  Ec},  {—Ec  <  E{ 0)  <  Ec}  or 
{P(0)  >  Ec}. 

The  local  relation  (6)  quantifies  the  hysteretic  rela¬ 
tion  between  E  and  P  for  homogeneous,  single  crystal 
compounds  having  uniform  effective  field  Ee  =  E.  To 
extend  this  mesoscopic  model  to  macroscopic  regimes 
involving  nonhomogeneous,  polycrystalline  compounds 
with  variable  effective  fields,  we  consider  certain  param¬ 
eters  to  be  manifestations  of  underlying  distributions 
rather  than  constant  values. 

Effective  field  effects  (see  [1])  are  incorporated  by  as¬ 
suming  that  the  field  at  the  lattice  level  is  normally 
distributed  about  the  applied  field  E  with  the  density 

f(Ee)  =  Cie-(^-s)2/5  (8) 

where  ci  and  b  are  positive  constants.  To  incorporate 
variations  in  the  lattice  structure,  we  assume  that  the 
local  coercive  field  Ec,  specified  by  (7),  is  either  nor¬ 
mally  or  lognormally  distributed.  In  the  former  case, 
it  will  exhibit  a  density  /  analogous  to  (8)  whereas  in 
the  latter  case,  it  has  the  density 

f(Ec)  =  C2e®MSc/.Ec)/2&]2 


E _ PrT) 

r]—2Yp^e  r\—2Yp^e 
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where,  if  b  is  small  compared  with  Ec,  Ec  denotes  a 
mean  coercive  field  at  which  dipoles  switch. 

The  macroscopic  polarization  model  is  then 


P(E)  =  C 


P(E,  e;  EC1  £)f(Ec)f{Ee)  dEedEc 


(9) 

where  P  is  given  by  (6).  It  is  detailed  in  [18]  that 
the  formulation  (9)  provides  an  energy  basis  for  certain 
Preisach  representations  with  the  difference  that  tem¬ 
perature  and  frequency  dependencies  are  incorporated 
in  the  kernel  of  (9)  rather  than  the  parameters  as  is  the 
case  for  Preisach  models. 

The  elastic  constitutive  relation,  in  the  absence  of 
internal  damping,  is  determined  from  the  equilibrium 
conditions 


d2G 

de2 


>  0 


which  yields 


cr  =  Ype  -  YpjP. 
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To  incorporate  Kelvin- Voigt  damping,  we  posit  that 
stress  is  proportional  to  a  linear  combination  of  strain, 
strain  rate  and  polarization  which  yields  the  1-D  con¬ 
stitutive  equation 

a  =  Ype  +  cpi  -  Yp1P  (10) 

where  cp  denotes  the  damping  parameter  at  fixed  po¬ 
larization.  In  combination  with  (9),  the  constitutive 
relation  (10)  quantifies  the  material  behavior  in  oper¬ 
ating  regimes  dominated  by  uniaxial  dynamics  as  is  the 
case  for  PZT  rods  employed  in  the  x-  and  y-stages  of 
several  present  AFM  designs. 

2.2.  2-D  Constitutive  Relations 

The  constitutive  behavior  of  a  PZT  shell  employed 
for  nanopositioning  differs  from  that  of  the  rod  in  two 
fundamental  aspects:  (i)  the  longitudinal  actuation  is 
due  to  g?3i  rather  than  (I33  mechanisms,  and  (ii)  lon¬ 
gitudinal  and  circumferential  stresses  and  strains  are 
coupled  due  to  the  curvature  (e.g.,  see  [2]).  The  na¬ 
ture  of  the  actuation  simply  yields  a  different  electrome¬ 
chanical  coupling  coefficient  which  we  denote  by  3  in 
this  case.  To  designate  the  coupled  material  behavior, 
we  let  ex ,  (7 x  and  eg ,  erg  respectively  denote  the  normal 
strains  and  stresses  in  the  longitudinal  and  circumferen¬ 
tial  directions  and  we  denote  shear  strains  and  stresses 
by  exg  and  axg.  Finally,  we  let  v  denote  the  Poisson 
ratio  for  the  material. 

To  simplify  the  discussion,  we  consider  the  develop¬ 
ment  of  constitutive  relations  in  the  absence  of  internal 
damping  and  refer  the  reader  to  [2,  7,  15]  for  the  exten¬ 
sions  necessary  to  incorporate  Kelvin- Voigt  damping. 
With  Yp  again  denoting  the  Young’s  modulus,  appro¬ 
priate  2-D  constitutive  equations  are 

Yp  Yp  3 

- n{ex  +  veg)  -  - — — P(E ) 

1  —  v*  1  —  1/ 

Yp  Yp3 

og  =  z - j  (e®  +  Vex )  ~  i - ~P(E) 

Yp 

G x6  —  @6x  —  777 1  \  \^x6 

2(1  +  v) 

/»oo  /»oo 

P{E)  =  C  /  P(E,e-,Ec,S)f(Ec)f(Ee)dEedEc 

J  0  J  —  00 

(n) 

with  P  specified  in  (6).  These  relations  are  employed 
when  constructing  the  shell  models  in  Section  3.2. 

3.  System  Models 

It  was  noted  in  Section  1  that  two  actuator  geome¬ 
tries  commonly  employed  for  nanopositioning  are  piezo¬ 
ceramic  rods  and  piezoceramic  shells.  In  this  section, 
we  construct  system  models  for  these  configurations 
based  on  the  1-D  and  2-D  constitutive  relations  devel¬ 
oped  in  Section  2. 


3.1.  Stacked  Actuator 


When  modeling  the  stacked  actuator  employed  in  the 
x-  and  y-stages  of  an  AFM,  we  make  the  assumption 
that  the  cross-sectional  area  A  is  small  compared  with 
the  length  i  and  we  consider  only  longitudinal  displace¬ 
ments  u  having  the  spatial  coordinate  x.  In  accordance 
with  present  stage  design,  one  end  of  the  rod  is  as¬ 
sumed  fixed  while  the  other  encounters  resistance  due 
to  the  connecting  mechanisms.  We  assume  that  this 
contribution  can  be  modeled  as  a  damped  elastic  sys¬ 
tem  with  mass  Ml,  stiffness  Aq.  and  damping  coefficient 
cl-  The  density  of  the  stacked  actuator  is  denoted  by 
p  and  the  Young’s  modulus  and  Kelvin- Voigt  damping 
coefficients  are  again  denoted  by  Yp  and  cp . 

Force  balancing  along  the  stacked  actuator  yields  the 
relation 


d2u  dJ\f 

pAW  =  77 


(12) 


where  the  resultant  J\f  =  JA  crdA  is  given  by 


A'='p4+^S-^p(ji» 

once  the  linear  relation  £  =  §y  is  employed  for  the 
strains.  The  nonlinear  and  hysteretic  map  between  in¬ 
put  fields  E  and  the  polarization  P  is  specified  by  (9). 
The  fixed-end  condition  yields  u(t,  0)  =  0  and  balanc¬ 
ing  forces  at  x  =  £  yields 

cP‘r\i 

A f(t,£)  =  - kLu{t,i )  -  d,  —  (t,£)  -  ML-^(t,t). 
Finally,  initial  conditions  are  taken  to  be  u(D,x)  = 

f(  0,7  =  0. 

The  relation  (12)  along  with  the  boundary  conditions 
provides  a  strong  form  of  the  stacked  actuator  model. 
For  numerical  implementation,  it  is  advantageous  to 
reduce  smoothness  requirements  on  approximating  el¬ 
ements  through  consideration  of  a  weak  form  of  the 
model  obtained  either  through  Hamiltonian  principles 
or  integration  by  parts.  Details  regarding  the  deriva¬ 
tion  of  a  weak  form  for  an  analogous  magnetostrictive 
rod,  and  a  description  of  the  resulting  ODE  system  ob¬ 
tained  through  a  finite  element  discretization  of  this 
weak  model,  are  provided  in  [5]. 

The  performance  of  the  model  is  illustrated  in  Fig¬ 
ure  4  where  the  model  prediction  is  compared  with  qua¬ 
sistatic  AFM  data  collected  at  0.28  Hz.  For  the  model 
construction,  normal  densities  /  and  /  were  employed 
when  quantifying  the  coercive  and  effective  field  dis¬ 
tributions.  It  is  observed  that  the  model  accurately 
quantifies  the  hysteresis  loss  and  hence  will  eliminate 
the  delay  associated  with  uncompensated  hysteresis  in 
present  control  designs. 


3.2.  Cylindrical  Actuators 

Secondly,  we  consider  the  development  of  a  model 
for  the  cylindrical  transducer  depicted  in  Figure  2(b). 
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Mx,Ms  and  Mxg  are  moment  resultants.  The  bound¬ 
ary  conditions  at  the  fixed  end  x  =  0  are  taken  to  be 
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Figure  4.  Quasistatic  experimental  data  from  the 
AFM  stage  depicted  in  Figure  2(a)  and  model  predic¬ 
tion. 


We  focus  on  the  actuator  employed  for  ^-displacements 
since  real-time  control  of  this  component  is  required 
to  maintain  constant  forces  between  the  sample  and 
micro-cantilever.  The  mass  of  the  shell  employed  for 
x-y  translation  is  combined  with  the  mass  of  the  sample 
to  provide  an  inertial  force  acting  on  the  free  end  of  the 
^-actuator. 

For  modeling  purposes,  we  assume  that  the  shell  has 
length  £,  thickness  h ,  and  radius  R.  The  axial  direction 
is  specified  along  the  ~-axis  and  the  longitudinal,  cir¬ 
cumferential  and  transverse  displacements  are  respec¬ 
tively  denoted  by  u,  v  and  w.  The  density  and  Young’s 
modulus  are  respectively  designated  by  p  and  Yp ,  and 
the  region  occupied  by  the  middle  surface  of  the  shell  is 
specified  by  To  =  [0,  l\  x  [0,  27t]  .  We  consider  the  case  in 
which  the  bottom  edge  of  the  shell  (a;  =  0)  is  clamped 
and  the  opposite  end  ( x  =  £)  is  acted  upon  only  by  the 
inertial  force  associated  with  the  combined  mass  m  of 
the  x-y  actuator  and  the  sample.  Finally,  we  simplify 
the  discussion  by  considering  the  model  development 
in  the  absence  of  internal  damping.  The  extensions 
necessary  to  incorporate  Kelvin- Voigt  damping  can  be 
incorporated  by  assuming  that  stress  is  a  linear  combi¬ 
nation  of  strain,  strain  rate  and  polarization  as  detailed 
in  [2,  15]. 

As  detailed  in  [2] ,  force  and  moment  balancing  yield 
the  Donnell-Mushtari  shell  equations 


R  hd2u  dAfx 

Rphdt2  R  dx 


dAfxg 
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Rph 
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~dP 


Rph 
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-  R 


dx2 


dAfxg  _ 
dx 

1  d2Me 
R  de 2 


2  AAxg 
dxdO 


+  Afg  =  0 


where  Afx,Afg  and  Afxg  are  general  force  resultants  and 


dw 

u  =  v  =  w  =  — —  =  0, 
dx 


and  the  conditions 
.  ,  d2U 

V.  =  -mw 

1  dMxs  _ 
Jx+  R  dd 


Afxt 


Mx 


R 


=  0 


Mx  =  0 


are  enforced  at  x=i.  The  first  resultant  condition  incor¬ 
porates  the  inertial  force  due  to  the  mass  of  the  piezoce¬ 
ramic  cylinder  employed  for  y-z  translation  along  with 
the  mass  of  the  sample. 

The  force  and  moment  resultants  are  specified  by  in¬ 
tegrating  the  stress  relations  (11),  or  the  product  of  the 
stress  and  moment  arm,  through  the  thickness  of  the 
shell.  In  the  absence  of  shear  stresses,  this  yields 

Afx  —  - - y(ex  +  veg)  — - - P(E) 

1  —  v2  1  —  v 


.r  yPh  ,  ,  Yph/3 

Afg  —  z - n(eg  +  vex )  — - P(E) 


Afxg  = 


1  —  v2 

Yph 

2(1  +  z/) 


l-i/ 


&x0 


and 


Mx  = 


Mg  = 


A4xg  = 


Yph3  .  ,  Yph3/3 
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vkx)  ~  yxyj - zP(E) 
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The  midsurface  strains  and  changes  in  curvature  are 
given  by 


ee  = 


&x6  — 


du 

d2w 

dx  ’ 

Kx  - 

dx2 

1  dv  w 

Kg  = 

1  d2w 

R  dd 

R  ’ 

R2  de 2 

dv 

1  du 

r  =  - 

2  d2w 

dx 

R  d6  ’ 

R  dxdd 

(13) 


As  with  the  stacked  actuator,  it  is  advantageous  to 
consider  a  weak  form  of  the  model  when  developing 
approximation  techniques.  Details  regarding  this  de¬ 
velopment  can  be  found  in  [2,  7]. 
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4.  Concluding  Remarks 

In  this  paper,  we  have  developed  1-D  and  2-D  consti¬ 
tutive  relations  quantifying  the  nonlinear  and  hystere¬ 
sis  behavior  of  PZT  transducers  employed  in  present 
nanopositioner  designs.  These  constitutive  relations 
are  then  combined  with  classical  rod  and  shell  theory  to 
provide  models  which  quantify  the  displacements  gen¬ 
erated  by  stacked  and  cylindrical  PZT  actuators.  This 
provides  a  design  tool  for  future  actuator  design  as  well 
as  a  framework  amenable  to  inversion  as  a  prelude  to 
control  design  employing  approximate  inverse  filters  to 
compensate  for  the  hysteresis  and  constitutive  nonlin¬ 
earities  inherent  to  the  materials  [12]. 
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